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We present a generalized Galitskii approach for the Bethe-Salpeter equation for the two-particle 
vertex function of a Fermi system with resonant interaction by accounting for the resonant state 
in the scattering potential and utilizing the universal form of the resonant scattering amplitude. 
The procedure can be carried out both for the normal as well as for the condensate state. In both 
cases, the vertex function in the vicinity of the resonance is shown to formally coincide with that 
obtained for a weakly attractive Fermi gas. Thus we justify the popular calculational framework in 
which results for the weakly attractive Fermi gas are formally extrapolated into the domain of strong 
coupling, and further to the repulsive side of the resonance, where molecular states are formed. 
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I. INTRODUCTION 

The problem of an interacting Fermi system with 
strong resonant coupling is relevant to many areas of 
physics, in particular, ultracold gasesi* 2 -' 3 - and high T c 
superconductors. 1 Experimental insight into this prob- 
lem has recently been obtained using ultracold atomic 
6 Li and 40 K gases ^ in which the s-wave scattering 
length a can be tuned from positive (repulsive) to nega- 
tive (attractive) values by crossing a Feshbach resonance, 
at which the scattering length diverges. These experi- 
ments started with the successful formation of long- living 
paired fermions in the strong-coupling regime,- which 
were subsequently seen to condense into a Bose-Einstein 
condensate (BEC) of bound molecular states.— Soon 
afterward the interaction was tuned to the attractive 
side^^ i 10 : 11 and it was established that the condensate 
crosses over to a BCS superfluid of extended Cooper 
pairs, 12 : 13 ! 14 ! 15 ' 16 : 17 

The first theoretical studies of the BCS-BEC crossover 
were initiated not long after the development of the BCS 
theory of superconductivity itself, when it was noticed 
that the nonzero solution of the BCS gap equation, com- 
plemented with the condition of a constant number of 
particles and regularized to eliminate an ultraviolet di- 
vergency, can be smoothly extrapolated to the domain 
of positive scattering lengths (without diverging at the 
resonance), where it describes molecular states! 18 : 19 ' 20 ! 21 

This mean-field interpolation between the two types of 
condensates is believed to be most accurate at T = 0, 
where bosonic excitations of paired fermions (Cooper 
pairs or molecular states of finite momentum) are not 
populated. Subsequent works established that the pair 
excitations induce fluctuation corrections which modify 
the thermodynamics especially in the strong-coupling 
limit ,22i23i24i25i26i27 ^ man y Q £ these works the fluctu- 
ation corrections are treated within the self-consistent 
T-matrix approximation which involves the ladder ap- 



proximation for the vertex function. 

The calculations beyond the mean field approxima- 
tion in the BCS theory were initiated by Gor'kov and 
Mclik-Barkhudarovj 2 ^ similar types of corrections were 
also applied to the BCS-BEC crossover.— Other re- 
cent approaches combine the ladder approximation with 
the fully self-consistent solution for the single particle 
Green function^ and account for the scattering of the 
weakly bound composite bosons. It should be noted 
that the scattering length for the two molecules in vacuo 
can be calculated exactly^ providing a useful guid- 
ance for the many-body calculations. Furthermore, the 
BEC-BCS crossover can also be studied numerically in 
Monte Carlo simulations ; 32 ! 33 ' 34 or via the renormaliza- 
tion group analysis with e and l/N expansions at the 
unitary point i 3 ^ 3 ^! 3 ^^ 8 .! 3 ^ 

Comparison with experiment and Monte-Carlo calcu- 
lations reveal that the ladder approximation captures 
the qualitative features for various physical quantities 
in the BCS-BEC crossover far better than it could be 
expected , 40 ' 41 This is also astonishing since the corre- 
sponding Bethe-Salpeter equation is conventionally de- 
rived by assuming the constant approximation for the 
scattering matrix^ which does not explicitly account 
for the resonant state. This approximation is strictly 
valid only for a weakly attracting Fermi gas. In order to 
study the strong-coupling limit, the resulting expressions 
are formally extrapolated, and these expressions are also 
used on the repulsive side of the resonance. Curiously, 
they then capture molecular states of the correct energy, 
and also recover the scattering length for the molecules 
in the Born approximation. Moreover, the results turn 
out to be nondivergent at resonance. These observations 
have motivated several works which address the justi- 
fication of solving the Bethe-Salpeter equation for the 
strong coupling case and in the presence of bound molec- 
ular states. In the context of high-T c two-dimensional 
superconductors the validity of the extrapolation pro- 
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cedure for the strong coupling regime was discussed by 
Randeria et a/4£, while for the ultracold Fermi gases the 
properties of the resonant scattering amplitude were used 
by Combescot^ to search for the universal form of the 
scattering matrix which would be applicable close to the 
resonance while recovering both limiting cases far from 
it. 

In this work we solve the Bethe-Salpeter equation for 
the two-particle vertex function of a Fermi system near 
the resonance, and show that the resulting expression 
coincides with the result of the extrapolation and regu- 
larization procedure. Our derivation proceeds via an ex- 
tension of the Galitskii formalism, 42 which reformulates 
the Bethe-Salpeter equation by using the scattering am- 
plitude instead of the scattering potential. In its original 
version, the Galitskii formalism is unsuitable to describe 
resonant scattering, since it does not explicitly account 
for the resonant state and as a result treats the case of 
positive scattering length a as a repulsive Fermi system. 
In our derivation we do not use any assumptions beyond 
the condition that the resonant state induces a broad 
Feshbach resonance. The state is eliminated by using 
the completeness and orthogonality of the set of all scat- 
tering states. This procedure can be carried out both for 
the normal state as well as for the condensate state, and 
leads to a fully analytical solution of the Bethe-Salpeter 
equation which is inherently valid in the vicinity of the 
resonance, as well as in the domain of the positive scat- 
tering length. The resulting expressions of this formal 
derivation turn out to be identical to the BCS extrapo- 
lation scheme with the renormalized contact interaction. 
Furthermore, our derivation constitutes a reorganization 
of the ladder approximation which demonstrates that the 
relevant expansion parameter does not diverge at reso- 
nance. Thus we establish a firm basis for the common 
conceptual framework which is used both in the mean 
field description of the BCS-BEC crossover, as well as in 
many calculations accounting for the fluctuation correc- 
tions. 

This paper is organized as follows. Section [II] defines 
the microscopic Fermi model on which all considerations 
arc based. In Section HTT1 we present the modified Galit- 
skii formalism for the Fermi gas in the normal state. In 
Section IIVI we present the derivation for the Fermi sys- 
tem in the condensate state, and also briefly discuss the 
gap equation. The results of this paper are summarized 
m Section El 



standard single-channel Hamiltonian 



2m 



M V'a(r) 



H = J2 [ dxUi(r) 

a J < 

+ \T.J drdr'i;i(r)i;l(r')U a p(r-r')^(r')i, a (r), 



(1) 

where ip a (v) is a fermionic field operator for a particle 
with spin a. In the following we work in the units m = 1, 
h = 1, and also set Boltzmann's constant ks = 1. The 
effective scattering potential U a p(r) describes resonant 
scattering of particles of opposite spin in the s-wave chan- 
nel. In the vicinity of a broad resonance with effective 
interaction radius ro <C kp 1 (where kF is the Fermi mo- 
mentum) , the scattering amplitude in the s-wave channel 
takes the universal form 4 ^ 



/(k) « - 



1 



77 + i|k| 



(2) 



where 77 = a -1 is the inverse scattering length. The scat- 
tering length a diverges at resonance, rj = 0. For 77 < 
the resonant level lies in the continuum, and the inter- 
action is attractive. For 77 > the resonant level turns 
into a bound state, and the interaction is repulsive. In 
the BCS limit 77 — > —00, the interaction mediates the 
formation of Cooper pairs, which are weakly bound in 
momentum space. In the BEC limit r\ — > +00, the inter- 
action mediates the formation of molecular states, which 
are weakly bound in real space. Because these notions 
refer to the condensate phase, we more broadly speak of 
the Fermi limit for 77 — * — 00, and of the Bose limit for 
j] — > +00. 

It should be stressed that the Hamiltonian (fl]) is 
purely fermionic, and can be contrasted to the more de- 
tailed Fermi-Bose models which explicitly account for a 
bound molecular state in the closed channeh 25 i 46 ' 47 i 48 
A great deal of work has been carried out in the 
past to give a detailed justification of the validity of 
both models and, in particular, to show their equiv- 
alence in the strong-coupling regime of a broad Fesh- 
bach resonance i 49 i 50 ' 51 i 52 ' 53 i 54 The universality arises for 
a combination of two facts. Firstly, at 77 = the scat- 
tering length diverges, which renders the Fermi energy 
Ep = kp/2 as the only relevant energy scale for the 
thermodynamics^ Secondly, in this unitary limit the 
contribution of the closed channel becomes negligible, 
which has been directly tested in experiments on 6 Li, 56 
even though the situation in 40 K may be a matter of 



II. SINGLE-CHANNEL MODEL FOR THE 
RESONANTLY INTERACTING FERMI GAS 



III. MODIFIED GALITSKII FORMALISM FOR 
RESONANT SCATTERING 



In order to describe an ultracold gas of fermionic atoms 
in the vicinity of a Feshbach broad resonance we use the 



We first assume that the Fermi gas is in the normal 
state. In the ladder approximation, the vertex function 
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T(pi,p2',P3,P4:) is then determined by the Bethe-Salpeter 
equation 



r(pl,P25P3,P4) =«(P1 -P3) / (^p"^ 

x G(pi - g)G(p 2 + q)T(pi - g,_p 2 + q;p 3 ,P4), (3) 

where spin indices are omitted for clarity, pi — (pi t o, p) 
is the 4D energy-momentum vector of an incoming (i = 
1,2) or outgoing (i = 3,4) particle, u(q) is the momen- 
tum representation of the interaction potential, and 



G(q) 



1 



?0 - £q ' 



Mi 



(4) 



is the Green function of a noninteracting Fermi system. 
We use the notation q = (qo, q), where q denotes the 
three-dimensional momentum vector and qo = iirT(2n + 
1) is the Matsubara frequency. 

In solving Eq. (J3J) one follows the conventional 
strategy^ 2 -, developed for the T = case, where the po- 
tential it(q) is eliminated in favor of the vacuo vertex 
function Tq, defined as the solution of Eq. © for two 
particles in absence of all the other particles, i.e., by set- 
ting fi — 0. An extension of this strategy for the finite 
temperature case is conveniently achieved by a subse- 
quent analytic continuation of the energy argument in 
T as 



E = g + 2 f i- g 2 /4. 



(5) 



Hereafter, we denote the total momentum of the scat- 
tered particles as g = pi + P2 = (ffo,g), while p — 
(pi — pa)/2 and p' — (p^ — Pa) 1 2 are the relative mo- 
menta before and after the scattering event, respectively. 

Eliminating it(p) in Eq. ([3]), the Bethe-Salpeter equa- 
tion takes the form 



r(p,p', 5 ) = r (p,p',s) 

d 3 q 



(2tt) 3 

The kernel K is defined as 



r (p,q, ff )^(q,.g)r(q,p', 5 ). (6) 



K(q,g) = tJ^G (I -q,flo- qo) G 



q,qo 



qo 
1 



E - 



(7) 



Using Eq. (fj| one obtains an explicit expression of the 
kernel, 



n F (£+) + n F (g-) 
E-a 2 



(8) 



where njr(e) = [exp(e/T) + 1] 1 is the Fermi function 
and 



5(f^) ; 



(9) 



In the Galitskii approach (see, e.g., Ref. |42|), the vacuo 
vertex function Tq is found by relating it to the two- 
momentum scattering amplitude /(p, k). However, the 
original formalism does not account for a possible bound 
molecular state in the scattering potential. In the follow- 
ing we assume that the scattering potential u(q) allows 
for a single bound state, with eigenfunction <p(x) and 
eigenvalue A. 

In the first step, Tq is written as an integral 



f d 3 q 

r (p,p',s) = J ^-^u(q)x(p- q,p,s), 



(10) 



where the scattering state x(PjP'jS) satisfies the equa- 
tion 

/d 3 
t^t3 »(q)x(p - q, p',g) 



(2n) 3 (E -p 2 + iO)J(p-p'). 



(11) 



The left-hand side of this equation is equivalent to a 
Schrodinger equation in momentum space. Therefore, 
x(PiP'>d) can b e written using the complete set of solu- 
tions of the Schrodinger equation, 



r d 6 } 

x (p,p',.g) = (£-p' 2 + *0) J 
E - p' 2 + iO 



d 3 k MpWUp') 



E-2\ + i0 



(2?r) 3 E - k 2 + iO 
P)^(p'). (I 2 ) 



where V'k(p) is the Fourier transform of the scatter- 
ing state ipk(r) with momentum k. The last term of 
this expression explicitly accounts for the bound state 
4>(p). The states V'k(p) can be expressed via the two- 
momentum scattering amplitude /(p,k), 



Mp) = (2tt) 3 5( P - k) + 



/(P,k) 



k 2 - p 2 + iO ' 
Inserting this expression into Eq. (|12[) one obtains 

C (7 3 k 

x(p,p',9) = <Mp) + J ^^k(p)r( P ',k) 



(13) 



1 



1 



+ 



E - k 2 + iO k 2 - p' 2 - iO 
E-p' 2 + iO 



E-2\ + i0 



(pW(p')- 



(14) 



After substituting Eq. (0 into Eq. flTUJ) one obtains Tq 
as a sum of two terms 

To (p, p', 9) = if (p, p', g) + T b (p, p', g) . (15) 

The first term in Eq. (fT5j) recovers the original Galitskii 
result 

/^ 3 k 
^33 /(P,k) f* (p'.k) 



(16) 



E - k 2 + iO k 2 - p' 2 - iO 
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The second term in Eq. (fT5|) originates from the bound 
state and is given by 



Il(p,p',3) = ?^ / T^uCqMp-q^V). 



E-2X J (2tt) 3 



(17) 



This expression still explicitly depends on the bound 
state, and at first glance it appears that detailed knowl- 
edge of <^>(p) is necessary in subsequent calculations. 
However, 4>(p) can be eliminated by relating it to /(k, p). 
To establish this relation we utilize the Schrodinger equa- 
tion for 0(p) in the momentum representation 



(2A - p 2 )0(p) 



/ ™«(q>(p - q) (18) 



and the completeness relation 

/r/ 3 k 
^3^k(p)V4(p') + HpW( P ') = (2nfS(p p'). 

(19) 

Multiplying both side of Eq. (jTTJ) by <f>(p), then using 
Eq. (fl"9|) and the definition (p~3j) one obtains the correc- 
tion of the Galitskii result due to the contribution of the 
molecular state as 



rg (p,p',9) = 



E - p' 2 2A - p 2 
E - 2A p 2 - p' 2 - 



it) 



d 3 k 

(2^)3 



/(P,k) 



x/* (P',k) 



k 2 -p' 2 -iO k 2 -p 2 + i0 

/(p,p')-r(p',p)V (20) 



Equations (|T5|) . (|T6|) and (|20j) reduce the calculation 
of the vacuo vertex to that for the two-momentum scat- 
tering amplitude /(p,k), defined by Eq. (fT3")) . This am- 
plitude is obtained by solving the Lippmann-Schwingcr 
equation. In the weak-scattering limit, /(p, k) can be ap- 
proximated by a constant 4tt / 77. Following Eq. (fT5|) . this 
leads to a constant approximation for To, as well. Since 
the vertex function then does not possess any poles, the 
constant approximation cannot be used to describe the 
bound states. In a more accurate analysis, applicable 
close to the resonance as well as in the limits r\ — ► ±00, 
one can utilize the fact that in the vicinity of a reso- 
nance the two-particle scattering amplitude for s-wave 
scattering /(p,k) assumes as similar universal form as 
the single-particle scattering amplitude /(k)^ 



/(p,k) 



4-7T 



ilkl 



(21) 



where in general r\ depends on p and k. For small mo- 
menta, r\ approaches a constant equal to the inverse scat- 
tering length. The leading corrections are quadratic, 
oc k 2 and cx p 2 , but the coefficients for these terms are 
proportional to the effective radius of the scattering po- 
tential tq i 42 ' 45 Under the condition of a broad resonance 



one has r^kp <C 1. Therefore, these terms can be ne- 
glected and 77 in Eq. (|2"Tj) can be assumed constant, equal 
to the inverse scattering length. On the molecular side of 
the resonance, where the potential permits a bound state, 
the inverse scattering length is furthermore related to the 
bound state energy via 2 A = —r\ 2 ^ 

It should be noted that Eq. (|21[) reduces to the cor- 
rect constant expressions in the weak-coupling limits r\ — > 
±00, and also satisfies the optical theorem; as we will see 
later, these observations are strongly linked to the suc- 
cess of the conventional cxtrapolation-and-regularization 
procedure for the vertex function. 

Substituting Eq. (J^TJ) into Eq. (JTSJ) when 77 < and 
into Eqs. (HHJ) and (f2"0"|) when r\ > and using 2 A = 
— 7y 2 one obtains a simple expression for the vacuo vertex 
function as 



T (E) 



An 



7] + iy/E' 



(22) 



where the branch cut of the square root lies on the real 
semi-axis E > 0, such that on the first Riemann sheet 
iy/—l = — 1. Recalling the definition of E one sees that 
the vacuo vertex function only depends on the total mo- 
mentum and energy of the molecular pair. Equation (|22p 
is applicable on both sides of the resonance and does not 
have any singularity at the resonance. On the repulsive 
side rj > 0, Eq. (|2"2"j) has a pole at E = —r\ 2 , which de- 
scribes a bound state with binding energy rj 2 . On the 
attractive side r\ < 0, the pole is absent: it is located on 
the unphysical second Riemann sheet of complex energy 
E. At resonance, Eq. (|22|) yields a square root singular- 
ity, which is different from the simple pole obtained in 
Ref. 0. 

For comparison it is instructive to recalculate Tq for 
a repulsive potential rj > using the original Galitskii 
result (fT6|) . i.e., without the bound state. In this case 
one obtains a different expression 



r (p',E) =4tt 



2t] 



1 



7] — iy~E 



(23) 



Contrary to Eq. (|22|) , Eq. (|23| does not have a pole that 
corresponds to the bound state. Equations (|2"2"|) and (|2"3"|) 
explicitly distinguish between potentials that have the 
same scattering length but do or do not permit a bound 
state, respectively. We note that the series expansion of 
Eqs. (f2"2")) and (f2"3")) with respect to large r\ differs only in 
the third-order term. Therefore, both Eqs. (|22"|) and (|23|) 
yield the same result when the Bethe-Salpeter equation 
<j6j) is solved in a second-order approximation, as done in 
the original Galitskii formalism^ 

We now proceed with the solution of Eq. (|6|) for the 
two- particle scattering amplitude (|22| . In doing so, we 
perform the analytical continuation defined by Eq. ([5]). 
The solution can then be written as 



T(g) = {T^+K(g)Y 



An 



77 + i\/E + AnK(g) ' 



(24) 
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where 



K(g) = 



d 3 q tanh(£ + /2T) + tanh(£_/2T) - 2 



(2tt) 3 



2(q2 - E) 



(25) 



In the limit ?/ — * oo, fi — > — oo, where the system is domi- 
nated by the molecular states, the vertex (j2~4")l reduces to 
the vacuo vertex function (f22|) . 

In the final step of the derivation we integrate the last 
term in Eq. (|25|) by parts. This cancels the term iy/E in 
Eq. ([M]), yielding 



r(s) 



47T 



7, + if (g) ' 



(26a) 



where 



f d 3 q J tanh(£ + /2T) + tanh(£_/2T) 2 



47T 2 \ 



(26b) 



This expression coincides precisely with the widely used 
extrapolated weak-coupling result, 1,2 including the ultra- 
violet regularization for the contact interaction (see, e.g., 
Ref. [H). The derivation procedure above demonstrates 
that in the ladder approximation, this expression remains 
strictly valid in the vicinity of the resonance, as well as 
on the molecular side 77 > 0. The reason for the fact that 
the vertex function (|26p . originally derived in the limit 
r/ — > — oo, correctly describes a molecular bound state 
for 77 > lies in the fact that T approaches Tq in Eq. 
(f2"2")l . There is also another, more intuitive reason why 
the formal extrapolation of the results for the weakly 
interacting Fermi limit to the domain 77 > correctly de- 
scribes the case with a molecular bound state and not the 
weakly repulsive Fermi gas. For a purely repulsive poten- 
tial without such a state, the resonant limit 77 — > can 
only be realized when the effective radius of the potential 
also becomes infinite. This would be in conflict with the 
continuity of the vertex function (|26[) as a function of 77, 
which also holds across the resonance at 77 = 0. 



IV. VERTEX FUNCTION IN THE 
CONDENSATE STATE 



The modified Galitskii approach described above can 
be equally applied to solve the Bethe-Salpeter equation 
for the vertex function in the condensate system. The 
condensate system has normal as well as anomalous ver- 
tex functions (with four incoming or outgoing lines) and 
the resulting Bethe-Salpeter equation becomes a system 



of equations 60 

/<i 3 q 
7^-TT u (q) 
.,„ 

{g(pi - q)G(p2 + q)T(pi - q,p 2 + q;P3,P4,) 

+ F( Pl - q)F(p 2 + q)Q( Pl - q,p 2 + q;p 3 ,p 4 )Y (27) 



d 3 q 



?! - q,P2 + q;p3,P4) 



{G(q-p 1 )G(-p2-q)e( Pl 
+ F*( Pl - q)F*{p 2 - q)T( Pl -q,p 2 + q;p 3 ,Pi)}, (28) 

The normal and anomalous Green functions G and F, 
respectively, are defined as 



G(q) 



90 +£q _ 

ql-^ * W qi-Al 



F(q) 



(29) 



where A is the BCS single particle spectral gap, related 
to the condensate density, and 



A q = + 



(30) 



while £ q is given in Eq. (f5|). Solution of this system 
along the lines of the Galitskii formalism outlined in the 
previous section yields the algebraic system 

T(g) = r„(ff){l - K(g)T(g) - S(g)Q(g)}, 

SO?) = -T (-g){K(-g)Q(g) + S*(g)T(g)}, (31) 

where the kernel K (g) is defined as the integral of Eq. Q 
over k and the Green function, G, defined by Eq. 
while S(g) is 



d 3 q 
(2tt) 3 



(32) 



The algebraic pair of Eqs. (|3Tj) is solved by 



r(s) = ^r-j , &{g) = „^ , (33) 

where the common denominator is given by 



2(.9) = (f? + WE + 4irK(g))(r] + iV E + AitK{-g)) 
- 16tt 2 S(g)S*(g). (34) 

As in the normal case the partial integration in the nor- 
mal kernel removes i\f~E and the result for the vertex be- 
comes equivalent to that in the weakly coupled neutral 
BCS system, which again by the virtue of the derivation 
is valid for arbitrary 77, including 77 = 0. 

We conclude the analysis by a brief discussion of the 
gap equation, which determines the order parameter A. 
According to the Thouless criterion, the stability of the 
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condensate requires that the bosonic excitation spectrum 
defined by the poles of T is gapless. This yields 



V 



E + AnK(0) + 4tt5(0) = 0. 



(35) 



[The second possible equation with the minus sign in 
front of S(0) yields instabilities in the corresponding 
bosonic excitation spectrum.] The explicit form of Eq. 



-V 



dq 



2A n 



■ tanh 



A, 
2T 



(36) 



This is the familiar BCS gap equation for a weakly inter- 
acting BCS system. However, as the underlying vertex 
function (f3"3"]l . Eq. (|3"6")l is valid (within the mean-field 
approximation) in the vicinity of the resonance as well 
as on its repulsive side, r\ > 0. The observation that at 
r\ > this mean-field gap equation, combined with the 
condition of the constant particle density, correctly de- 
scribes the bound molecular states (see, e.g., Refs.flaand 
l2lh is again ensured by the correct limit T — ► Tq of the 
vertex function in the molecular limit rj — > oo, fi — ► — oo. 



V. CONCLUSIONS 

In this work we have presented a procedure of solv- 
ing the Bethe-Salpeter equation for the two-particle ver- 
tex functions in the ladder approximation in the vicin- 
ity of a broad resonance, as encountered for ultracold 
Fermi gases close to a Feshbach resonance. In order to 
do this we have extended the Galitskii formalism to ac- 
count for the molecular states in the interaction potential 
and also utilized the universal form of the resonant scat- 
tering amplitude in the vicinity of the resonance. This 
allowed for the exact solution of the Bethe-Salpeter equa- 
tion at resonance and its vicinity. It is also valid in the 
weak-scattering limit, including the limit of the molecular 
states. The solution is regular in the BCS-BEC crossover 
and coincides exactly with the standard results obtained 
for the weakly attractive Fermi gas. This agreement pro- 
vides a link between the resonance regime and the weak- 
interaction regime and serves to justify the widely used 



calculational framework in which the weak-interaction re- 
sults (with ultraviolet regularization) are extrapolated to 
the BCS-BEC crossover and further to the domain of 
molecular states. 

The derivation of the vertex function can be carried out 
both in the normal as well as in the condensate state. As 
the latter is used to extract the equation for the conden- 
sate density (the gap equation) , the extrapolation proce- 
dure for the BCS gap equation onto the molecular side of 
the resonance is also perfectly justified (within the lad- 
der approximation). Thus our results establish the ro- 
bustness of the phcnomcnological extrapolation from the 
BCS limit used in many previous works. 

We have not discussed the validity of the ladder ap- 
proximation itself, which is obviously questionable when 
the system is close to the resonance. The applicability of 
the ladder approximation is conventionally established 
for a small gas parameter kp/\r}\ <C 1. This condition 
apparently fails at resonance, where the gas parameter 
diverges. Still, it has been found that the ladder approxi- 
mation combined with a self-consistency procedure deliv- 
ers reliable qualitative, and to some extent even quanti- 
tative information in the complete crossover regime (see, 
e.g., Refs. Q] and@). The presented solution of the Bethe- 
Salpeter equations in the vicinity of the resonance reveals 
that the relevant small parameter in the resonant system 
is cx kp\f(kp)\, where is the absolute value of 

the scattering amplitude. This quantity coincides with 
the conventional gas parameter in the weak-interaction 
limit, and, although not small, is not divergent at res- 
onance, approaching a constant of order unity instead. 
While not providing a small expansion parameter, this 
observation provides further insight into the apparent 
successes of the ladder approximation in studies of the 
resonant scattering systems. 
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